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Abstract 

The tensor virial theorem for subsystems is formulated for three- 
component systems and further effort is devoted to a special case 
where the inner subsystems and the central region of the outer one 
are homogeneous, the last surrounded by an isothermal homeoid. The 
virial equations are explicitly written under the additional restrictions: 
(i) similar and similarly placed inner subsystems, and (ii) spherical 
outer subsystem. An application is made to hole + vortex + bulge 
systems, in the limit of flattened inner subsystems, which implies three 
virial equations in three unknowns. Using the Faber-Jackson rela- 
tion, R e oc Oq, the standard Mh-cto form (Mh oc Oq) is deduced from 
qualitative considerations. The projected bulge velocity dispersion to 
projected vortex velocity ratio, r] = (o- B )33/{[(«v)gg] 2 + [(o"v) gq ] 2 } 1/2 , 
as a function of the fractional radius, yev = Rb/Rv, and the frac- 
tional masses, tobh = Mb /Ma and mev = Mb /My, is studied in 
the range of interest, < myH = Mv/Mh < 5 (Escala, 2006) and 
229 < tobh < 795 (Marconi and Hunt, 2003), consistent with obser- 
vations. The related curves appear to be similar to Maxwell velocity 
distributions, which implies a fixed value of rj below the maximum 
corresponds to two different configurations: a compact bulge on the 
left of the maximum, and an extended bulge on the right. All curves 
lie very close one to the other on the left of the maximum, and parallel 
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one to the other on the right. On the other hand, fixed msH or tobv, 
and 2/BV; are found to imply more massive bulges passing from bot- 
tom to top along a vertical line on the (Ot/bv 7 ?) plane, and vice versa. 
The model is applied to NGC 4374 and NGC 4486, taking the frac- 
tional mass, 7«bh, and the fractional radius, ubv, as unknowns, and 
the bulge mass is inferred and compared with results from different 
methods. In presence of a massive vortex (myH = 5), the hole mass 
has to be reduced by a factor 2-3 with respect to the case of a massless 
vortex, to get the fit. Finally, the assumptions of homogeneous inner 
bulge and isotropic stress tensor are discussed. 

keywords - Black hole physics; galaxies: bulges. PACS codes - 98. 62. J s 



1 Introduction 

The existence of large masses confined in a restricted central region of galaxies 
(hereafter quoted as "holes" ) , is widely supported by current high- resolution 
observations (for a review see e.g., Ferrarese and Ford, 2005; Merritt, 2006). 
There is an amount of increasing evidence, that compact objects at the centre 
of (sufficiently large) galaxies are supermassive black holes (e.g., Maoz, 1998; 
Miller, 2006). In addition, the presence of self-gravitating nuclear disks (with 
typically several hundreds parsecs in radius; hereafter quoted as "outer vor- 
texes") seems necessary to continuously replenish inner Keplerian accretion 
disks (hereafter quoted as "inner vortexes"), allowing substantial hole mass 
growth (e.g., Escala, 2006). Hole growth can successfully be modelled as a 
rapid succession of vortexes with comparable mass, related to a merger event 
(King et al., 2008). Most of the gas residing in the outer, self-gravitating 
disk, is expected to be either turned into stars, or expelled by the first gen- 
eration of high-luminosity stars, on a rapid (almost dynamical) time-scale 
(King et al., 2008). 

In the following, the disk made of the inner vortex and the outer (self- 
gravitating) vortex, shall be quoted as the "vortex". Masses up to about 
one fifth the inner bulge (i.e. bulge stars inside the isophote with equato- 
rial axes coinciding with the edge of the vortex) have been deduced from 
observations (e.g., Downes and Solomon, 1998). The occurrence of vortexes 
is expected to be a common feature during the evolution of galaxies host- 
ing holes. Accordingly, central regions of large galaxies may be modelled as 
the superposition of three main components: a central hole, a surrounding 
vortex, and an embedding bulge. 

In this view, the tensor virial theorem for subsystems (Brosche et al., 
1983; Caimmi et al., 1984; Caimmi and Secco, 1992; Caimmi, 2007) could 
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be a useful tool of investigation. In the case under discussion there are 
nine virial equations, three for each subsystem. In absence of a satisfactory 
theory of quantum gravity, the mass distribution within the hole cannot be 
described, which reduces the number of virial equations from nine to six. 
The further restriction to an axisymmetric gravitational potential, implying 
an axisymmetric matter distribution for each subsystem, further reduces the 
number of virial equations from six to four. The final restriction to flattened 
vortex and hole configurations, reduces the number of virial equations from 
four to three. 

Keeping in mind that the virial theorem deals with global instead of local 
properties of matter distributions, simple even if unrealistic density profiles 
can be used for investigating main features and general trends. To this aim, 
the hole and the vortex are modelled as concentric and coaxial, flattened, 
homogeneous spheroids, the inner bulge as a homogeneous sphere, and the 
outer bulge as an isothermal spherical corona, which allows an analytical 
formulation (Caimmi and Secco, 2002). 

To close the system of three equations, the number of unknowns must be 
reduced to three. The variables are: three mass-averaged square velocities, 
three masses, and three equatorial semiaxes. Hole mass and equatorial semi- 
axis can be related via the expression of the gravitational (or Schwartzschild) 
radius (e.g., Landau and Lifchitz, 1966, Chap. XI, §97). Typical values of 
bulge and vortex equatorial semiaxes, can be deduced from observations (e.g., 
Downes and Solomon, 1998). The vortex to hole mass ratio cannot exceed 
a value of about 5 (e.g., Escala, 2006), deduced from observations (Downes 
and Solomon, 1998). An empirical correlation seems to exist between bulge 
and hole mass (e.g., Marconi and Hunt, 2003). Mass-averaged bulge square 
velocities may safely be approximated by square luminosity-weighted second 
moments of the line-of-sight velocity distribution within the half-light radius, 
which is weakly dependent on the details of the orbital distribution (Cap- 
pellari et al., 2006). In fact, a similar quantity is also used for studying the 
correlation between hole mass and bulge velocity dispersion (e.g., Gebhardt 
et al., 2000). Taking the above informations into due account, the number 
of unknowns in the virial equations may be reduced from nine to three. 

In addition, the combination of two virial equations yields a relation be- 
tween mass-averaged velocity ratios, mass ratios, and equatorial axis ratios, 
regardless of the values of mass-averaged velocities, masses, and equatorial 
semiaxes. For assigned bulge to vortex and vortex to hole mass ratios, and 
vortex to hole equatorial axis ratio, the mass-averaged bulge to vortex ve- 
locity ratio is a function of the bulge to vortex equatorial axis ratio, and the 
related trend can be studied. In particular, the bulge mass can be deter- 
mined by the knowledge of other parameters, and the related value can be 
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compared with results from different methods. 

The current paper is organized as follows. The virial theorem for subsys- 
tems, with regard to three-component systems, is formulated in Sect.EJ The 
model and its application to hole + vortex + bulge systems are the subject 
of Sect. El where the results are also presented and discussed, with regard to 
the special cases of NGC 4374 (M84) and NGC 4486 (M87). The conclusion 
is drawn in Sect.Hl 

2 Basic theory 

A general theory of two-component matter distributions has been exhaus- 
tively treated in earlier papers (e.g., Caimmi and Secco, 1992; Caimmi and 
Marmo, 2003; Caimmi and Valentinuzzi, 2008) and an interested reader is 
addressed therein and in parent investigations (MacMillan, 1930, Chap. Ill, 
§76; Limber, 1959; Neutsch, 1979; Brosche et al, 1983; Caimmi et al, 1984). 
The related tensor virial theorem is expressed by 3 x 2 independent equa- 
tions, each containing 2 + (2 — 1) terms corresponding to the kinetic energy 
and the self potential energy of the subsystem under consideration, and to 
the tidal potential energy due to the gravitational potential induced by the 
other subsystem, respectively. 

In dealing with iV subsystems, the tensor virial theorem is expressed by 
3xN independent equations, each containing 2+(iV — 1) terms corresponding 
to the kinetic energy and the self potential energy of the subsystem under 
consideration, and to the tidal potential energies due to the gravitational 
potential induced by the other subsystems, respectively. Accordingly, the 
results found for two-component systems can be used, by inserting in the 
virial equations N — 2 additional terms, related to the tidal potential energies 
of the extra components. 

To this aim, the formulation of earlier attempts (Caimmi and Secco, 1992, 
2002) shall be used. If not otherwise stated, for sake of convenience, mat- 
ter distributions shall be conceived as continuous media instead of discrete 
particle sets (e.g., Limber, 1959; Caimmi, 2007). In the following, general 
definitions and related explicit expressions will be presented. Readers mainly 
interested to an application of astrophysical interest (hole + vortex + bulge 
systems) might directly go to the next Section [31 

2.1 The tensor virial theorem for subsystems 

Let i, j, k, denote the subsystems of a three-component matter distribution. 
Owing to the additivity of the gravitational potential, the potential-energy 
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tensor and the potential energy may be cast into the form: 

(la) 

U U V 

E po t = ^2{E U ) SC \ + X ~~ <5«u)(-^™)int ! (lb) 
p = l,2,3 ; g = l,2,3 ; u = i,j,k ; uv = ij,ik,jk,ji,ki,kj ; (lc) 

where <5 U „ is the Kronecker symbol, [(£ , u )sei]pg an d (-Eu)sei the self potential- 
energy tensor and the self potential energy, respectively, [(E u ) int \ pq and (E u ) int 
the interaction potential-energy tensor and the interaction potential energy, 
respectively. The related explicit expressions are (e.g., Caimmi and Secco, 
1992): 

[(E u ) scl ] pq = / pu(x 1 ,x 2 ,x 3 )x p — -d 3 S u 

= ~\ I Pu(x 1 ,x 2 ,x 3 )[V u ( y x 1 ,x 2 ,x 3 )} pq d 3 S u ; (2a) 
(E u ) scl = I p u (xi,x 2 ,x 3 )J2x s ^-d 3 S u 

J S u s=l (J^s 

1 r „ 
= ~o Pu(x 1 ,x 2 ,x 3 )V u (x 1 ,x 2 ,x 3 )d S u ; (2b) 

£ J S u 

[{E uv ) int \ pq = ~ ( p u (x 1 ,x 2 ,x 3 )[V v (x 1 ,x 2 ,x 3 )] pq d 3 S u ; (3a) 

^ J S u 

(E uv ) int = ~f p u {x 1 ,x 2 ,x 3 )V v {x 1 ,x 2 ,x 3 )d 3 S u ; (3b) 
^ J s u 

where S is the volume, p the density, V pq and V the gravitational tensor 
potential and potential, respectively (Chandrasekhar, 1969, Chap. 2, §10). 

The tensor and the scalar virial theorem for a single subsystem within 
the tidal field induced by the other ones, read: 

2[(-E'u)kin]pg + [(E u ) se \\p q + [(E uv ) t id] pq + [(E uw ) tid ] pq = ; (4a) 

2(-E«)kin + (-E'u)sei + (E uv ) tid + (E uw ) tid = ; (4b) 
u = i,j,k ; v = j,k,i ; w = k,i,j ; (4c) 

where [(-Eu)kmW an d are the kinetic-energy tensor and the kinetic 

energy, respectively, [{E uv ) t i d ] pq and (-E^tid the tidal potential-energy tensor 
and the tidal potential energy, respectively. The related explicit expressions 
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are (e.g., Caimmi and Valentinuzzi, 2008 

-I 

2Js, 



[(E u )un] Pq = 77 / Pu{xi,x 2 ,x 3 ){v u ) p {v u ) q d 3 S u ; (5a) 



If 3 

(E u ) kin = - p u (x 1 ,x 2 ,x 3 )J2[( v u)s] 2 d 3 S u ; 



(5b) 



[(E uv ) tid ) pq = Pu(xi,x 2 ,x 3 )xp-— !L d 3 S u ; (6a) 

J S u 0Xq 

(E uv )tid= p u (x 1 ,x 2 ,x 3 )J2 x s-^- 1( i 3 S u ; (6b) 



where (v u ) p (v u ) q = (v u ) p (v u ) q and [(f u ) s ] 2 = [(f u ) s ] 2 are arithmetic means cal- 
culated within the infinitesimal volume element, d 3 ^ = dxi dx 2 dx 3 , placed 
at P(xi,x 2 ,x 3 ) (Binney and Tremaine, 1987, Chap. 4, §4.3). The tidal poten- 
tial energy may be conceived as the virial of the u component in connection 
with the tidal field induced by the v component (Brosche et al., 1983). 

The tensor and the scalar virial theorem, expressed by Eqs. (I4al) and (14bl) . 
are the extension of earlier results related to one-component (e.g., Chan- 
drasekhar, 1969, Chap. II, §11; Binney and Tremaine, 1987, Chap. 4, §3) and 
two-component (e.g., Caimmi et al., 1984; Caimmi and Secco, 1992) systems. 

The tensor and the scalar virial theorem for the whole system read: 

2(-E'kin)pq + (E pot ) pq = ; (7a) 
2E kin + E pot = ; (7b) 

which can be made explicit using Eqs. (llal) - ffTbl) . An alternative formula- 
tion can be obtained summing Eqs. (|4lil) - (l4b|) over all the subsystems. The 
combination of the related expressions yields: 

tid]p5 [(-^mi)int]p<? 

(E )tid + (E vu ) ti d — (E uv ) int + (E vu ) int ; (8b) 
uv = ij, ik,jk ; vu = ji, ki, kj ; (8c) 

according to earlier results related to two-component systems (Caimmi and 
Secco, 1992). 

In the general case of subsystems with different shapes, the tensor and 
scalar virial equations, Eqs. fj4"al) and (|4bl) . need numerical integrations to be 
made explicit, except for concentric and coaxial homogeneous ellipsoids, one 
completely lying within the other (e.g., Caimmi and Secco, 1992). The ef- 
fect of different shapes is larger for flat density profiles, with respect to their 
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monotonically decreasing counterparts. Explicit formulae also exist if the 
outer ellipsoid, in turn, is surrounded by a similar and similarly placed het- 
erogeneous homeoid, where the inner surface coincides with the boundary of 
the outer ellipsoid (Caimmi and Secco, 2002). In the more realistic situation 
of homeoidally striated ellipsoids (Roberts, 1962), a further restriction aimed 
to simplify the calculations, is that different subsystems are also similar and 
similarly placed (e.g., Caimmi, 1993; Caimmi and Marmo, 2003; Caimmi and 
Valentinuzzi, 2008). 

In dealing with an application of the tensor virial theorem to a three- 
component system of astrophysical interest, the choice of the related density 
profiles shall be limited to the above mentioned special cases. 



3 The model 



3.1 The hole, the vortex, and the bulge 

Aiming to a simple application of the tensor virial theorem, inner regions of 
large galaxies shall be idealized as three-component systems, one completely 
lying within the other: the hole, the vortex, and the bulge, respectively. The 
outer vortex is self-gravitating, and the presence of the inner vortex allows 
mass accretion inside the last stable orbit. The inner bulge is enclosed by 
the isopycnic (i.e. constant density) surface with equatorial axes coinciding 
with the edge of the vortex. 

In absence of a satisfactory theory of quantum gravity, the mass distri- 
bution within the hole cannot be described, which reduces the number of 
tensor virial equations, Eqs. (I4al) . from nine to six. The further restriction to 
an axisymmetric gravitational potential, implying an axisymmetric density 
profile, for each subsystem, further reduces the number of tensor virial equa- 
tions from six to four. The final restriction to flattened vortex configurations, 
reduces the number of tensor virial equations from four to three. 

The hole shall be modelled by a flattened homogeneous spheroid potential, 
where the equatorial axis coincides with the radius of the particle inner stable 
circular orbit of the related Schwartzschild black hole: 

Rr = 3i? s ; R g = — ^ — ; (9) 

where G is the constant of gravitation, c is the light velocity in (baryonic) 
vacuum, M is the mass, R g is the gravitational radius (e.g., Landau and Lif- 
chitz, 1966, Chap. XI, §97), and the index, H, denotes the hole. The choice of 
a Plummer potential (Escala, 2006) would imply a more complicated formu- 
lation (Caimmi and Valentinuzzi, 2008). The selected value for Rn, expressed 
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by Eq. ([H]) and related to Schwartzschild's metrics, makes an acceptable 
compromise between a corotating particle inner stable circular orbit radius, 
i?H = (l/2)-R s , and its counterrotating counterpart, Ru = (9/2)R g , related 
to Kerr's metrics. For the photon inner stable circular orbit, Ru = (3/2)R g in 
Schwartzschild's metrics, and Ru = (l/2)R g (corotating), i?H = 2R g (coun- 
terrotating), in Kerr's metrics. For further details refer to specialized reviews 
(e.g., Bilic, 2006). 

The vortex may be conceived, to a first extent, as self-gravitating even 
if it is still not clear if fragmentation or turbulence are dominating (e.g., 
Escala, 2006). For simplicity, both the vortex circular velocity, v, and the 
inner bulge velocity dispersion, a, are assumed constant with radius, which 
is only strictly valid for an isothermal sphere (Escala, 2006). Even within the 
most massive gaseous vortexes known so far, the gravitational potential is 
dominated by the inner bulge. With regard to the inner bulge, the following 
values are deduced from observations (Downes and Solomon, 1998): 

Mi , \ 

™ IN = — ^5; ( 10 ) 

where My and M\ are the mass of the gaseous vortex and stellar inner bulge, 
respectively. 

The circular velocity at the edge of the vortex, is obtained by the condition 
of centrifugal equilibrium, as: 

2 Cv GMt Ch Ghh GM, 
= ™ + ™ P~ + Cl ~R~ ; (U) 

777.IN -Kv ^IH -Kv -n-V 

where ( is a shape factor and, in particular, ( = 1 for spherical-symmetric 
mass distributions and ( = tt/2 for flattened homogeneous spheroids. 
The inner bulge velocity dispersion at the edge of the vortex reads: 

( a i) = ; ( 12 ) 

provided the density profile declines as in an isothermal sphere. 
The combination of Eqs. fllOI) - f|T2l yields: 



K) 2 = ^ + — + (^i) 2 « (^i) 2 ; (is) 




which holds a fortiori for less massive vortexes, miN > 5. 

The validity of the vortex steady state model (the so called a disk) as a 
zeroth order approximation, has successfully been tested by recent numerical 
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experiments (Escala, 2006). The vortex shall be represented as a flattened 
homogeneous spheroid potential, similar and similarly placed with respect to 
the hole. 

Strictly speaking, the mass inside the hole volume should be subctracted 
which, in the case under discussion, amounts to: 

AM V = (^) 3 M v ; (14) 



(15) 



\Ry/ 

and the ratio of subtracted to remaining mass reads 

AM V (Rn/Ryf 



Mv-AMy 1-ORhARv) 3 ' 
which shows that the subtracted mass, AMy, may safely be neglected. 

The inner bulge, < r < Ry, shall be modelled by a homogeneous sphere 
potential, and similar considerations can be made for the subctracted mass: 

AM 1 =(^f) 3 M 1 ; (16) 



Ri 

keeping in mind that R\ = Ry. 

The outer bulge, Ry < r < Rb, shall be modelled by a truncated isother- 
mal sphere potential. The related density profile (including the inner bulge) 
is: 

Pb = Pv/b(£b) ; £b = -jr- ; (z/Bvr 1 = ^ ; (17a) 

Kb -Kb 

f it \ _ j 1 ; < £ B < (z/bv)" 1 ; / 17K x 

/B(ee) - j d/Bv6)- 2 ; (2/bv)- 1 < £b < 1 ; (17b) 

where, in general, tnb = Ry and pnb = Pv are a scaling radius and a scaling 
density, respectively. 

In conclusion, the central part of a large galaxy shall be modelled as an 
inner, flattened homogeneous spheroid (the hole), embedded within a similar 
and similarly placed, flattened homogeneous spheroid (the vortex) which, in 
turn, is lying within a homogeneous sphere with equal equatorial semiaxes 
(the inner bulge), surrounded by a concentric, isothermal, spherical corona 
(the outer bulge). 



3.2 The tensor virial equations 

The four tensor virial equations which, among the nine expressed by Eqs. (T4]), 
are relevant to the application of interest, are: 

2 [(£v)kin] pp + [(£v)sei] pp + [(£vhU pp + [(£v B )tid] pp = ; (18a) 

2 [(£ B )kin] pp + [(£ B )sel] pp + [(£BH)tid] pp + [(£ B v)tid] pp = J (18b) 
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where 23 has been chosen as symmetry axis, which implies p — 1,3, and all 
the terms may be written explicitly, under the assumptions previously made. 
The vortex and the bulge kinetic-energy tensors read: 



2[(B 
2[(E : 



V kin 



PP 



B ikin 



Vy « CT B 



PP 





Mpp 


2 

+ 


( a v) pp 






2 

+ 





V 



1,3 ; 



(19a) 

(19b) 
(19c) 



where 5 is the Kronecker symbol, v pp and (7 pp £11*6 the circular velocity and 
the random velocity components along the corresponding axis, respectively. 

In the special case of homeoidally striated ellipsoids (Roberts, 1962), the 
self-potential energy tensors read (e.g., Caimmi, 1993): 



sol 



PP 



-{vx) sel —^{tix)p , 



R 



x 



B. n 



^pq 



£p2£p3 



:i + s)- 3/2 (i + er ! .s i 



c 2 

-pl c 



X = V,B ; (20a) 
-^(l + dar^da ; (20b) 



P 



1,2,3 
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1,2,3 



p^i^r ; (20c) 



where i? is the major equatorial semiaxis, (^x) se i are profile factors, B p shape 
factors which, for axisymmetric configurations, may be analytically expressed 
(e.g., Chandrasekhar, 1969, Chap. 3, §21; Caimmi, 1991, 1995), and e pq axis 
ratios. For further details refer to AppendixlAl 

The vortex self potential energy profile factor reads (Caimmi and Secco, 
1992): 

K) sel = ^ ; (21) 

and the shape factors are expressed in AppendixlAl 

The bulge self potential energy profile factor reads (Caimmi and Secco, 
2002)0: 

3 15(?/bv) 2 - 1%bv - 103/BV ln(y B y) . ^ 



v^B, 



scl 



10 



(3j/bv - 2) 2 

related to the density profile, expressed by Eq. (JT7J). 

In the general case of triaxial configurations, the tidal potential-energy 
tensors depend on the density profiles and on the shape of the ellipsoids. 
The isopycnic surface of the outer ellipsoid which is still embedding the inner 



lr The related formula in the parent paper, Eq. (51) therein, is different due to printing 
errors. 
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one, is tangent to the latter boundary at a top axis, (a»)t. For oblate-like 
configurations, the inner more flattened than the outer, (a i ) i = (a$)i. 

The mass of the outer ellipsoid surrounded by the isothermal homeoid, is 
(Caimmi and Secco, 20020 

M B = M mas (M B ) ; (23a) 

(M B ) = — (pB)o(aB)i(aB)2(aB)3 ; (23b) 

N mas = (!/Bvr 3 [l + 3(!/BV-l)] ; (23c) 

Vbv = (jtev)t = 5 ( 23d ) 

where (pb)o is the central density, (Mb)o the mass of a homogeneous ellipsoid 
with same boundary and density equal to the central density, and the major 
axis of the inner bulge and the vortex coincide. Accordingly, Eq. ( 123ajl takes 
the equivalent form: 

M B = M I [l + 3(y B v-l)] ; (24) 

where Mi is the mass of the inner bulge. 

In the special case of two concentric and coaxial homogeneous ellipsoids, 
the more flattened one completely lying within the other with a tangential 
point at the top major equatorial semiaxis, Rx, surrounded by a isother- 
mal (X=V) or homogeneous + isothermal (X=H) homeoid where the inner 
surface coincides with the boundary of the outer ellipsoid, the tidal potential- 
energy tensors read (Caimmi and Secco, 2001, 2002): 



[(£xyWU = S Pq ^ M m jr^K(B Y ) p ; (25a) 

-ftx [{yYx)p\ 

GM 2 ( 5 1 

[(E Y x)m] Pq = -<W^f Mm { 2 [(2/yx)i] 2 Fy(2/yx) + ($yx) p } (B Y ) p (25h) 

I \ 3 m YX 1 M Y . , (a Y ) P ^ 1 fo - v 

^xYj tid = ^77 t^t — v — ; m Yx = ~~jT7~ > U/yx p = t — r > 1 ; 25c 

10 K2/yx)i] 3 W mas ^x {ax)p 

(ex) Pq _ (yyx) q . X = H,V; Y = B ; p=l,2,3; g = 1, 2, 3 ; (25d) 



(eY)p 9 (yyx) P 

^b I/bv = 2 7 ; ^bU/bh = t V5-7 ^ + 2 7 ^ ; 25e 

(Z/BV) 2 (?/Bv) 2 (Z/BH) 2 (l/Bv) 2 



2 With regard to the profile factor, the related formula in the parent paper, Eq. (49) 

£| has to be replaced by £| 



therein, is different due to printing errors. In addition, £| has to be replaced by £| in 



Eq. (50) therein. 
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1,2,3 



1,2,3 



;(25f) 

(25g) 
(25h) 



where z/ ti( j are profile factors and B pr are shape factors which, for axisym- 
metric configurations, may be analytically expressed (e.g., Chandrasekhar, 
1969, Chap. 3, §21; Caimmi, 1995). For further details refer to AppendixEl 
In the special case of two concentric and coaxial homogeneous ellipsoids, 
one completely lying within the other (e.g., Caimmi and Secco, 1992), the 
tidal potential-energy tensors read: 



[(£hv) 



tid 



[(3 



VHitid 



P<1 



pq 



GMl , , [WJ 



pi p 
GM? 



[By) 



p i 



l^HVJtid 

Fy(yvn) 



m VH 1 



- (^v) ti d [Mwa) + ($vh) p ] (B v ) p ; 



10 [(y 



(_y_ 

(y 



\ iui] 3 (vv)i 

1 2 



i ; 



VHJl 



VhJp_ 



VHjpJ 



(y 



1 2 



VHjp 



VHji 



(26a) 

(26b) 
(26c) 

(26d) 



which may be conceived as a three-component system, as discussed above, 
where the surrounding homeoid is homogeneous instead of isothermal or ho- 
mogeneous + isothermal. 

In the special case under discussion, the shape factors take the expression: 



($yx)i = 

($yx)3 = 

($vh)i = 
ch = e v 
( q v)p _ 
(a v ) r 



YXj2 
9 



1 + — (1 

10 v 



($vh) 2 = ($vh)s 
= e 5 e B = 1 ; 
(°h)p 



[an) 



(y 



VR)p 



(y 



VH r 



yvn 



(27a) 

(27b) 

(27c) 
(27d) 

(27e) 



and the four tensor virial equations of interest, Eqs. ( |T8i) . due to (|T9 l -( l27l) . 
may be written under the explicit form: 



[{vy) qq f + [(a y ) 



3 GM Y \ a 



10 Ry 



1 + 



1 



^vh(?/ 



VHj 



;(?/vh) s 



12 



2 m-Bv 



33J 



3 3y B v - 2 

£7 



l2 3 GMy 
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q = l,2; (28c) 



2 1 GM B / 15(?/bv) 2 - 1%bv - lOysv 1h(!/bv) 




where, in the special case of flattened configurations, e = 0, the shape factors 
reduce to a/e = n/2 and = 0, see AppendixlAl which yields (av)33 = 0, as 
expected. The validity of Eqs. (I28cj) and (128dl) implies an anisotropic bulge 
stress tensor, [(t>B) g(? ] 2 + [( cr B)gg] 2 < [(c"b)33] 2 , which is due to the presence 
of more flattened, inner subsystems. In fact, an anisotropic stress tensor is 
necessary to maintain a spherical shape in the case under discussion. 

The tensor virial equations, Eqs. (I28al) -( l28dl) . make a system of four equa- 
tions in thirteen unknowns: two rotation velocity components, (ity)ii = 
(vy) 2 2 and (vb)u = (^3)22; f° ur peculiar velocity components, (<xv)n = 
(<r v )22, (o"v)33, (cb)h = (0^)22, and (cr B ) 3 3; three masses, M H , My, and 
Mb; three major equatorial semiaxes, Rn, Ry, and Rb', and one axis ratio, 
e = en = ev- The hole mass, Mh, and the model hole major semiaxis, Ru, 
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are connected via Eq. which reduces the number of unknowns to twelve. 
Further reduction can be made using observational constraints and/or addi- 
tional assumptions. 

Rotation and peculiar velocity components could be deduced from obser- 
vations, together with the vortex and bulge major semiaxes. The remaining 
four unknowns would be three masses and one axis ratio. In the flat limit 
(e — > 0) Eq. (128bj) reduces to an indeterminate form, = 0, leaving three 
equations in three unknowns, Mh, My, and Mb. Of course, any other alter- 
native might be exploited, yielding a system of four (or three) equations in 
four (or three) unknowns. 



3.3 Application 



According to an earlier model of massive vortexes (Escala, 2006), the hole 
and the vortex mass shall be taken as Mh = lO 7 m and My = 5 • lO 7 m , 
respectively. A flat disk (e B = ey = 0) makes an acceptable approximation 
to an assumed major equatorial semiaxis, Ry = 125pc, and thickness, Zy = 
4pc, (Escala, 2006) implying ey = Zy/Ry = 0.032. The particularization of 
Eq. (E]) to the selected hole mass yields: i?H = 2.87 ■ 10~ 6 pc. 
The related parameters appearing in Eqs. (I28al) -( l28dl) are: 



myn 



y V H = 4.35-10' 







(29) 



and the substitution into the above mentioned equations, keeping in mind 
that Eq. (128bl) cannot be used for flat configurations, yields: 
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5 ln(j/Bv) + 



(30c) 



where the terms containing negative powers of ?/vh ^> 1 and y BB > yyn, have 
been neglected. 
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In the current model, the vortex is described as a homogeneous, self- 
gravitating, flattened spheroid, which implies the relations: 



1 



(E- 



Vjkin 



M 



V 



-M v (i?v) 2 fi v 



-M Y [v v (R 



,„, ; (31a) 
(31b) 



where Qy is the (solid-body) angular velocity, [(t>v) 2 ] the mean square ro- 
tation velocity, and Vy(Ry) = {2.5[(vy) 2 ]} 1 ^ 2 the rotation velocity at the 
edge. Strictly speaking, the presence of the hole implies differential rotation 
towards Keplerian velocities, but this effect can be neglected beyond the hole 
sphere of influence (e.g., Ferrarese and Ford, 2005). For the selected reference 
case, Eqs. fl3D) yield \(vy) qq ] 2 /[(a v ) 33 ] 2 « 0.4 or {[vy(Ry)) qq } 2 /[K) 33 ] 2 « 1, 
which is consistent with Eq. fflBl . as expected for massive inner bulges (e.g., 
Escala, 2006). 

In the limit of null vortex and hole mass, Eq. (I30cp reduces to: 



1 GM B 15(y BV ) 2 - Uy BV - 10y BV \n(y By ) 



5 R B (32/bv-2) 2 
and the combination of Eqs. (I30cl) and (1321 yields: 



(32) 



1 M H R b [(a B ; 
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m BB 
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C(m Bn ,y BV ) 



GM B 
3y BV - 2 



(o"bo) 33 



:>>:', 



({m B H,y B v) 



Vbv 



1 19 -i — 1 

- + (1 + m VH ) ln(yBv) + ^ m VH 



(33a) 
(33b) 



where in the range of interest, 1 < y B y < 10 4 , < my B < 5, the con- 
tribution of the factor, C( m BH,2/Bv), scales from about 0.05 to 2, related to 
(wivh, Vbv) — (5, 10 4 ), (0, 1), respectively. A selected reference case, (m V H, 2/bv) 
(5, 14), yielding M B = lO lo m for an assumed M B = 10 7 m Q , corresponds to 
C(5, 14) = 0.156690. In the limit of an infinitely extended bulge, y B y — > +oo, 
C( m BH,?/Bv) 0, due to Mb — > +oo, for the isothermal sphere. Then the 
fractional mass, m BB , to a first extent, may be considered as independent of 
the vortex mass and bulge to vortex major semiaxis ratio. 

The mere existence of a fundamental plane (Djorgovski and Davis, 1987; 
Dressier et al., 1987) indicates that structural properties in spheroids (el- 
liptical galaxies and spiral bulges) span a narrow range, suggesting some 
self-regulating mechanism must be at work during formation and evolution. 
In particular, projected light profiles exhibit a large degree of homogeneity 
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and may well be fitted by the r 1//4 de Vaucouleurs law. Accordingly, a nar- 
row range may safely be expected also for the factor, 1 — [(030)33/ (o"b)33] 2 5 
appearing in Eq. ( l33l) . which makes the hole mass, Mh, depend only on the 
product, Rb [(c"b)33] 2 , to a first extent. 

If, in addition, the bulge density profile is close to an isothermal sphere, 
then the central velocity dispersion projected on the line of sight, a , is close 
to the rms velocity dispersion averaged on the mass, ((73)33. Then the hole 
mass depends only on the product, R e &o, to a first extent. Using the Faber- 
Jackson relation, R e oc <Tq, the standard Mh-cto form (Mh oc Gq) is obtained 
(Escala, 2006). 

A different analysis can be performed fixing the bulge mass, Mb, or 
the bulge to vortex mass ratio, rriBv, and considering the velocity ratio, 
V — ( cr B)33/{[(^v)gg] 2 + [{ (T v)qq\ 2 } 1 ^ 2 1 as a function of the bulge to vortex 
major axis ratio, ?/bv- The related explicit expression is obtained by the 
combination of Eqs. ( 1271) . (I28aj) . ( 128dl) . particularized to flat inner subsys- 
tems, which represents isofractional mass (i.e. fixed wbv, ^bh, m vn) curves 
in the (Ch/BV??) plane. 

The fractional masses, by definition, are related as msv — ^bh/^vh, 
where < myH < 5 in the cases of interest. On the other hand, a cor- 
relation between bulge and hole mass, yielding 229 < ttibh < 795, with 
a preferred value, 771bh ~ 427, has been deduced from recent observations 
(Marconi and Hunt, 2003). The assumption of constant ?tt,bh implies isofrac- 
tional mass curves in the (Ot/bv 7 ?) plane depend on a single fractional mass, 
mvn. Different cases in the range of interest, for wbh = 427, are shown by 
full curves plotted in Fig. CD Counterparts corresponding to an estimated 
upper (msH = 229) and lower (ttibh = 795) limit to the Mh-Mb rela- 
tion (Marconi and Hunt, 2003), are shown by dashed and dotted curves, 
respectively. Different curves of each class (from top to bottom) corre- 
spond to myH = 0,1,2,3,4,5, respectively. The reference configuration, 
(wbh, ?tivh, i/bv) = (427,5, 14), is represented by an asterisk. 

An inspection of Fig. [TJ discloses the following features. 

(i) All curves exhibit an extremum point of maximum. Lower ttj-bh and/or 

myH correspond to lower maximum coordinates, and vice versa. 

(ii) All curves are close, one with respect to the other, when they are rising, 

on the left of the maximum. 

(iii) All curves are parallel, one with respect to the other, when they are 

declining, on the right of the maximum. 

(iv) Curves related to different fractional masses, wbh and myH, could be 

nearly coincident for a convenient choice of values. This is the case for 
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0.5 i i i i 

10 20 30 40 50 

Figure 1: Isofractional mass curves in the (Ot/bv 7 ?) plane for ttibh — 427 
(full curves), which represents the preferred value of a recent interpolation 
of the empirical M^-M-q relation (Marconi and Hunt, 2003). Counterparts 
corresponding to an estimated upper (wbh = 229, dashed curves) and lower 
(^bh = 795, dotted curves) limit to the Mh-Mb relation (Marconi and Hunt, 
2003), are also shown. Different curves of each class (from top to bottom) 
correspond to m V H = 0, 1, 2, 3, 4, 5, respectively. The reference configuration, 
( m BH) m VH;2/Bv) = (427,200, 14), is represented by an asterisk. 



(mBH,mvH) = (427, 1) and (795,4), plotted in Fig. CD as full and dotted 
lines, respectively. 

3.4 Special cases 

Aiming to see the model at work, two special cases shall be considered, 
where vortex rotation curves have been deduced from observations: NGC 
4374 (M84) and NGC 4486 (M87). Parameter values of interest are listed 
in Tab.[TJ Bulge data and related uncertainties are taken from a recent 
investigation (Cappellari et al., 2006), and errors are propagated using a 
standard quadratic formula. Vortex data are taken from different sources 
for NGC 4374 (Bower et al, 1998) and NGC 4486 (Macchetto et al., 1997). 
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NGC/M 



4374/84 



4486/87 



-R e /kpc 

a e /(kms _1 ) 

Ry/pc 

Vy sini/ (kms -1 ) 



6.15=Fl.88 (1) 

278^14 (1) 

164^50 (2) 

250^150 (2) 

80^5 (2) 

0.15SS (2) 



7.96^3.23 (1) 

298^15 (1) 

35^14 (3) 

500^100 (3) 

51^5 (3) 

0.32^0.09 (3) 



i/deg 
M H /M 10 



Table 1: Parameter values related to NGC 4374 (M84) and NGC 4486 (M87). 
Captions: R e - bulge effective (half-light) radius; a e - luminosity-weighted 
second moment of the line-of-sight velocity dispersion within the effective 
radius; R Y - vortex radius at the end of the rotation velocity profile; Vy 
- vortex rotation velocity weighted on the mass, MyVy = j s pvv^ ot dSy] 
i - inclination angle; Mr - hole mass; Mio = 10 10 m Q . See text for de- 
tails. The velocity conversion factors from [cmgs] to [kpc M 10 Gyr] system 
of measure and vice versa are: 1 kpc Gyr _1 =0.978461 942 118 946 6 kms -1 
and lkms -1 =1.022012 156 992443 kpcGyr" 1 , respectively References: (1) 
Cappellari et al. (2006); (2) Bower et al. (1998); (3) Macchetto et al. (1997). 

Rotation velocities weighted on the mass are roughly estimated as the mean of 
the maximum and the plateau value of the rotation curve, with no correction 
for inclination effects, whose contribution is thought to be negligible in the 
framework of the above approximation. The related uncertainty is assumed 
to be equal to the deviation from the mean. Radii correspond to the end of 
measured rotation curves. The related uncertainty is assumed to be equal 
(in percentage) to its bulge counterpart, <tr/ / R'y = <tr b /Rb- 



The virial equations, Eqs. fl30al) and fl30cl) . provide the remaining two 
unknowns: the fractional mass, rriBn, and the fractional radius, y^y. The 
resulting position of NGC 4374 and NGC 4486 on the (Ot/bv??) plane, is 
(Vbv,v) = (37.5,1.11) and (113.68,0.60), respectively. The dimensions of 
the error box can be evaluated by use of standard quadratic propagation 
formulae; the result is (2<7 yBN , 2a v ) = (32.34,1.44) and (130,0.30), respec- 
tively. The comparison with model predictions is shown in Figs.[2]and[3]for 
NGC 4374 and NGC 4486, respectively. The position of each galaxy at the 
centre of the error box is marked by an asterisk. The value of the fractional 
mass, wbh, can be read near the corresponding curve on the right of the 
error box. The value of the fractional mass, myn, is null in all cases except 
lower full, dashed, and dotted curves, where m V H = 5 and m BH remains 
unchanged with respect to their upper counterparts (value labelled therein). 
The extreme cases still compatible with the error box, are represented by 
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dot-dashed curves. The centre of the error box is fitted by wbh = 230 (NGC 
4371) and mBH = 153 (NGC 4486); the related curves are not plotted in 
Figs. [2] and El to avoid confusion. 

3.5 Discussion 

The application of the tensor virial theorem for subsystems to hole + vortex 
+ bulge systems implies a number of restrictive assumptions to gain intrinsic 
simplicity related to analytical formulation. This is why the hole, the vortex, 
and the inner bulge are modelled as homogeneous, and the outer bulge as 
homeoidally striated. Accordingly, the bulge mass, Mb, the inner bulge 
mass, Mi = M B (Ry), and the fractional radius, t/bv — Rb/Rv, are related 
byEq.flSD. 

For assigned fractional masses, ttibh and msv, the velocity ratio, r] = 
(o"B)33/{[("yv) w ] 2 + [(°v)gg] 2 } , versus the fractional radius, y B v = Rq/Rv, 
is represented by a selected curve on the (Ot/bv 7 ?) plane, as shown in Fig.[H 
Interestingly, a fixed value of 77 below the maximum corresponds to two 
different configurations: a compact bulge on the left of the maximum, and 
an extended bulge on the right. 

On the other hand, fixed mev and jjbv, imply larger wbh passing from 
bottom to top along a vertical line, t/bv = const, as shown by the intersections 
of the above mentioned vertical line with the corresponding curves, see e.g., 
2/bv = 50 in Fig.[TJ Accordingly, larger wbh imply larger values of rj and vice 
versa. Similar results hold for the fractional mass, m-Qv, keeping wbh fixed. 

Upper curves of each series are related to massless vortexes, myH = 0. 
Test particles moving on circular and coplanar orbits, within a few hundred 
parsecs from the central hole, may be considered as massless vortexes. Then 
the knowledge of orbital parameters together with bulge effective radius and 
velocity dispersion, defines the position of the system on the (0?/bv ? 7) plane. 
The isofractional mass curve where the representative point lies, in turn, 
defines the value of the fractional mass, ttibh- 

The application to NGC 4374 and NGC 4486, shown in Figs. and El 
yields 8.5 < m B H < 2250 and 33 < m B H < 400, respectively. The empirical 
M B - M H correlation (Marconi and Hunt, 2003), 229 < m BH < 795, lies 
within the uncertainty range for NGC 4374 but only partially for NGC 4486. 
The knowledge of the vortex density profile would provide a better determi- 
nation of the mass-weighted rotation velocity, (Vy) qq , and then reduce the 
error box. In presence of a massive vortex (myH — 5), the hole mass has to 
be reduced by a factor 2-3 to get the fit. 

The bulge mass can be estimated from Figs. [2] and El using hole mass 
values deduced from gas dynamics listed in Tab.[TJ The result is: Mb/Miq = 
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34.35+^:65 and 48.96+?^ for NGC 4374 and NGC 4486, respectively. 

On the other hand, the bulge mass may be inferred from the knowledge 
of the mass to luminosity ratio, the /-band magnitude, and the distance 
modulus (e.g., Caimmi and Valentinuzzi, 2008), using recent data and related 
uncertainties (Cappellari et al., 2006). The result is: Mb/Mi = 72.70=f67.50 
and 91.94 =f 80.92 for NGC 4374 and NGC 4486, respectively. 

An alternative approach consists in dealing with the structural parameter 
estimate (Bender et al., 1992) of M B = (5/G)alR e , where a Q is the bulge 
central velocity dispersion projected on the line of sight, and R e is the bulge 
effective radius. An application to NGC 4374 can be found in an earlier 
attempt (Bower et al., 1998). Using recent data and related uncertainties 
(Cappellari et al, 2006) yields M B /M 10 = 72.70^22.62 and 108.12^47.84 for 
NGC 4374 and NGC 4486, respectively. The uncertainty should be multiplied 
by ( T , if 5(~ < 5 < 5( + is the scatter of the above relation. 

Finally, the substitution of hole mass values listed in Tab. [I] into the 
empirical M B -M H correlation, 229 < m BH < 795, with a preferred value, 
m BH = 427, (Marconi and Hunt, 2003), yields M B /M 10 = 64.05+?^ and 

136.64^SK for NGC 4374 and NGC 4486 > respectively. 

A comparison between the above results is shown in Fig.Hl Even if the 
current method is affected by a large upper uncertainty for NGC 4374, still 
it appears to be consistent with the other ones. 

In general, a selected point on the (Oi/bv 7 /) plane is crossed by infinite 
curves of the kind plotted in Fig.[H each characterized by a proper choice 
of fractional masses, (ttt, B v, tjbh) ■ 111 the limit of a massless vortex, m B y — > 
+oo, the model may be improved by replacing the inner bulge radius, R\ = 
Ry, with Ri = i?n- Accordingly, the bulge resembles a truncated isothermal 
sphere, where the central divergence of the density profile is avoided by the 
presence of a very small, homogeneous inner bulge. 

The bulge velocity dispersion, ((73)33, as a function of the bulge effective 
radius, R e , in the limit of a massless vortex, is plotted in Fig.[5]for Ri = Ry 
(full curves) and Rj = Ru (dashed curves), related to: (m B H,-Rv/pc) = 
(427, 125), (230, 164), (153, 70), in connection with Figs. d El El respectively. 
The position of the reference configuration (Fig.[T]) is marked by an asterisk. 
The position of NGC 4374 and NGC 4486 and the related uncertainties, taken 
or deduced from recent observations (Cappellari et al., 2006), are marked by 
asterisks and rectangles, respectively. It can be seen that the assumption of 
a homogeneous inner bulge makes bulge velocity dispersion systematically 
decrease for fixed bulge effective radius, but the error is comparable with 
observational uncertainties. 

The assumption of isotropic bulge stress tensor is acceptable for the spe- 
cial cases to which the model has been applied. More specifically, both NGC 
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4374 and NGC 4486 appear to be only slightly flattened. Even if the shape is 
entirely due to velocity dispersion, the stress tensor is expected to be moder- 
ately anisotropic provided the system is seen edge-on, and can be considered 
isotropic to a first extent. It is the case for NGC 4486, but not for NGC 
4374 (Bender et al., 1992). If, on the other hand, the system is seen head- 
on, the line of sight coincides with the polar axis and what is deduced from 
observations is related to ((73)33, regardless of the properties of the stress 
tensor. Strictly speaking, the degree of anisotropy and the inclination angle 
should be specified, in order to express the bulge velocity dispersion compo- 
nent along the polar axis, ((73)33, m terms of the luminosity- weighted second 
moment of the line-of-sight velocity distribution within the bulge effective 
radius, R e . 

4 Conclusion 

The tensor virial theorem for subsystems has been formulated for three- 
component systems and further effort has been devoted to a special case 
where the inner subsystems and the central region of the outer one are homo- 
geneous, the last surrounded by an isothermal homeoid. The virial equations 
have explicitly been written under the additional restrictions: (i) similar and 
similarly placed inner subsystems, and (ii) spherical outer subsystem. An 
application has been made to hole + vortex + bulge systems, in the limit 
of flattened inner subsystems, which implies three virial equations in three 
unknowns. 

Using the Faber- Jackson relation, (i?B) e oc ctq, the standard M H -(7o form 
(M H oc <7q) has been deduced from qualitative considerations. The projected 
bulge velocity dispersion to projected vortex velocity ratio, 77 = (cr B ) 33 /{[(vy) qq \ 
[(°"v)<?<?] 2 } 1//2 > as a function of the fractional radius, y^y = Rb/Rv, and the 
fractional masses, ttibh = Mb/ Mb and ttibv = Mb/ My, has been stud- 
ied in the range of interest, < m V H = M V /M H < 5 (Escala, 2006) and 
229 < wbh < 795 (Marconi and Hunt, 2003), consistent with observations. 

The related curves have been shown to be similar to Maxwell velocity dis- 
tributions, which implies a fixed value of 77 below the maximum corresponds 
to two different configurations: a compact bulge on the left of the maximum, 
and an extended bulge on the right. All curves have been seen to lie very 
close one to the other on the left of the maximum, and parallel one to the 
other on the right. 

On the other hand, fixed 777.33 or rriBy, an d 7/bv> have been found to imply 
more massive bulges passing from bottom to top along a vertical line on the 
(OyBvv) plane, and vice versa. The model has been applied to NGC 4374 
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and NGC 4486, taking the fractional mass, meH, and the fractional radius, 
tjbv, as unknowns, and the bulge mass has been inferred and compared with 
results from different methods. In presence of a massive vortex (myH = 5), 
the fit has been provided by hole mass 2-3 times lower with respect to the 
case of a massless vortex. Finally, it has been shown that the assumptions 
of homogeneous inner bulge and isotropic stress tensor hold to an acceptable 
extent, at least in the special cases taken into consideration. 
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Appendix 

A Ellipsoid shape factors: special cases 

Let di, a 2 , a 3 , a x > a 2 > a 3 , be semiaxes of a generic ellipsoid, and e pq = a p /a q 
related axis ratios. The shape factors, B p and B pr , defined by Eqs. (I20bl) and 
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Table 2: Values of ellipsoid shape factors, B p and B pr , related to limiting 
configurations defined by the values of the axis ratios, e 2 i and e 3 i. For ob- 
long configurations, Bi and Bu diverge as e 2 7, Byi and -B13 as e 2 . For flat 
configurations, -B13 and B 23 diverge as e -1 . 



( |25gP , respectively, in terms of their usual counterparts (e.g., Chandrasekhar, 
1969, Chap. 3, §21; Caimmi, 1991, 1995), read: 

Br, = e^e^Ar, ; (34) 

(35) 
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(36) 

j p -ft-p 

and a number of special cases are listed in Tab. [2] for limiting configurations. 

The range of validity of the related independent variables, €f, e p , e Q , and 
e, is: 



< e f < 1 ; 
e = e 3 i, < e < 1 



< e p < 1 



< e < 1 ; 



oblate 



.-1 



(37a) 
(37b) 

e = e 2 { = £31 , 1 < e < +00, prolate ; (37c) 
and the explicit expression of the functions, a, 7, wx, w 2 , w 3 , reads: 
e^[arcsiWl^_ e ] . oblate; 

(38) 

prolate; 
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where, in addition (e.g., Caimmi, 1991): 



w 3 = - (l - e 2 w 2 ) 



2a + 7 = 2 ; lim a = lim 7 = - ; lima = ; lini7 = 2 ; 

e— >1 e— >1 3 e ~ ^ e^O 



a 7T 

hm — = — ; 

6^0 e 2 



lim a = 1 ; lim 7 = ; lim £7 = 

e— >+oo 



lim e 7 = +00 ; lim w 2 = lim w 2 = - 

e— >1~ 5 



25 



2.0 



1.5 



1.0 



0.5 



0.0 

20 40 60 

y B v 

Figure 2: Comparison between model predictions and data from observations 
for NGC 4374 (M84). The position of the galaxy at the centre of the error box 
(deduced from values listed in Tab.[T]) is marked by an asterisk. The value 
of the bulge to hole mass ratio, mBH, can be read near the corresponding 
curve on the right of the box error. The value of the vortex to hole mass 
ratio, mvH, is null in all cases except lower full, dashed, and dotted curves, 
where m V H = 5 and ttt-bh remains unchanged with respect to their upper 
counterparts (value labelled therein). The extreme cases still compatible 
with the error box, are represented by dot-dashed curves. The centre of the 
error box is fitted by tjibyl = 230 (not represented). 
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Figure 3: Comparison between model predictions and data from observations 
for NGC 4486 (M87). The position of the galaxy at the centre of the error box 
(deduced from values listed in Tab.[T]) is marked by an asterisk. The value 
of the bulge to hole mass ratio, mBH, can be read near the corresponding 
curve on the right of the box error. The value of the vortex to hole mass 
ratio, mvH, is null in all cases except lower full, dashed, and dotted curves, 
where m V H = 5 and ttt-bh remains unchanged with respect to their upper 
counterparts (value labelled therein). The extreme cases still compatible 
with the error box, are represented by dot-dashed curves. The centre of the 
error box is fitted by tjibyl = 153 (not represented). 



27 



5 1 ' ' | 

4 |— * 1 

3 | * 

2 V — *— H 

I L x NGC 4374 (M84) J 

[ , , I 

200 400 600 

M E /M 10 



5 1 ' ' | 

4 1 — * i 

3 | x 

2 V * 

• L x NGC 4486 (M87) J 

I , , I 

200 400 600 

M B /M 10 



Figure 4: The bulge mass deduced from different methods, for NGC 4374 
(M84, top) and NGC 4486 (M87, bottom). Different values in ordinate are 
deduced from: 1 - the empirical Mb-Mh correlation, 229 < wbh < 795, 
with a preferred value, wbh = 427, (Marconi and Hunt, 2003), using hole 
mass values listed in Tab.[TJ 2 - the structural parameter estimate (Bender 
et al., 1992) of M B = (^/G)a^R e , where a a is the bulge central velocity 
dispersion projected on the line of sight, and R e is the bulge effective radius 
(the upper and lower limit should be multiplied by ( + and (~ , respectively, 
where < 5 < 5( + is the scatter of the above relation); 3 - the knowledge 
of the mass to luminosity ratio, the J-band magnitude, and the distance 
modulus (e.g., Caimmi and Valentinuzzi, 2008), using recent data and related 
uncertainties (Cappellari et al., 2006); 4 - Figs. [2] and El using hole mass 
values deduced from gas dynamics listed in Tab.[TJ Values of Mb/Mi for 
NGC 4374: 1 - 64.05+^1; 2 - 72.70^22.62; 3 - 72.70^67.50; 4 - 34.35^1;|. 
Values of M B /M 10 for NGC 4486: 1 - 136.64+?^;^; 2 - 108.12 t 47.84; 3 - 
91.94 =f 80.92; 4 - 48.96l?f£;jg. 
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Figure 5: The bulge velocity dispersion, (<7b)33, as a function of the bulge 
effective radius, R e , in the limit of a massless vortex, for R\ = Ry (full 
curves) and R\ = i?n (dashed curves), related to: (?tibh, -Rv/pc) = (427, 125), 
(230, 164), (153, 70), in connection with Figs.Cfl [21 [3X respectively. The po- 
sition of the reference configuration (Fig.[T]) is marked by an asterisk. The 
position of NGC 4374 (Fig.[2D and NGC 4486 (Fig.ED and the related uncer- 
tainties, taken or deduced from recent observations (Cappellari et al., 2006), 
are marked by asterisks and rectangles, respectively. 
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